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ABSTRACT 

Both microlensing surveys and radio-frequency observations of gas flow imply that 
the inner Milky Way is completely dominated by baryons, contrary to the predictions 
of standard cold dark matter (CDM) cosmology. We investigate the predictions of 
the Modified Newtonian Dynamics (MOND) formula for the Galaxy given the mea- 
sured baryon distribution. Satisfactory fits to the observationally determined terminal- 
velocity curve are obtained for different choices of the MOND's interpolating function 
fi(x). However, with simple analytical forms of fi(x), the local circular speed v c (Rq) 
can be as large as 220kms _1 only for values of the parameter ao that are excluded 
by observations of NGC 3198. Only a numerically specified interpolating function can 
produce v c (Rq) = 220 km s -1 , which is therefore an upper limit in MOND, while the 
asymptotic velocity is predicted to be v c (oo) — 170 ± 5 kms -1 . The data are probably 
not consistent with the functional form of fi(x) that has been explored as a toy model 
in the framework of Bekenstein's covariant theory of gravity. 

Key words: gravitation - Galaxy: kinematics and dynamics 



1 INTRODUCTION 



Recent advances in the modelling of the Milky Way indi- 
cate that, contrary t o the predictions of cold dark mat- 
ter (CDM) cosmology JPiemand. Moore fc Stadel20oJ. and 
references therein), the inner Galaxy is completely domi- 
nated by baryons. We ask whether current Galaxy models 
are more compatible with the Modified Newtonian Dynam- 
ics (MOND). 

The evidence for accelerating cosmic expansion can 
be interpreted either as signalling the presence of a large 
amount of dark energy (quintessence), or as an indication 
that Einstein's minimal theory of gravity must be modi- 
fied by adding a cosmological constant. Similarly, the evi- 
dence for flat galactic rotation curves is conventionally in- 
terpreted as evidence for massiv e halos o f cold dark mat- 
ter (CDM), but twenty years ago iMilgroml dl98otl suggested 
that flat rotation curves might signal the need to mod- 
ify Newtonian dynamics in regions where the acceleration 
is smaller than a critical value ao- This proposal, dubbed 
MOND, was then refined by the introduction of a non- 
relativistic field equation for the mod ified gravitational po- 
tential dBekenstein fc Milgromlll984h . There is now a sig- 
nificant body of evidence that, for whatever reason, there 
is a cha racteristic acceleration ap ~ cHo associated with 
galaxies JSanders fc McGaughlfeoofl . Furthermore, MOND 
has a remarkable ability to account for features in the phe- 
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nomenology of galaxies that were unknown when Milgrom 
introduced the theory. 

Recent developments in the theory of gravity have 
added plausibility to the case for modificatio n of gravity 
rather than addition of exotic matter. First, Bekenstein! 

has presented a Lorentz-covariant theory of gravity 
that has a MONDian behaviour in t he appropriate limit. 
Second, it has become recognized (e.g. Gripaios 2004) that 
spontaneous symmetry breaking in an effective field the- 
ory of gravity might well lead to loss of Lorentz invariance 
of the type required by Bekenstein. A great deal of work 
needs to be done to determine whether a theory such as 
that of Bekenstein is compatible with observations of the 
cosmic microwave b ackground and large-scale structure (e.g. 
ISkordis et a!l2 005l. But knowledge that it is possible to em- 
bed MOND within a physically acceptable dynamical frame- 
work, and that this framework is not unattractive from the 
point of view of mathematical physics, must make us take 
more seriously the phenomenological successes of MOND. 

On the other hand, there is growing evidence that the 
dark halos of galaxies do not conform to CDM simula- 
tions: they are cuspy in neither low-surface-brightness (e.g. 
Bosma 2004) nor high-surface-brightness galaxies, the cen- 
ters of which are completely baryon-dominated (Gerhard et 
al. 2001, Gentile et al. 2004, Cappellari et al. 2005). More- 
over the shapes of the rotation curves are intimately con- 
nected to the underlying stellar-light profile, just as modified 
gravity predicts. In this paper, we study the particular case 
of our own Milky Way galaxy, which is the best example of 
a baryon-dominated high-surface-brightness galaxy. 
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In Section 2 we review the evidence that the inner 
Galaxy is dominated by baryons. In Section 3 we inves- 
tigate the predictions of MOND for the Galactic circular- 
speed curve and the vertical equilibrium of the Solar neigh- 
bourhood. Section 4 sums up. 



2 THE BASEL MODEL OF THE MILKY WAY 

iBissantz fc Gerhard! i2002h (hereafter BG) present a model 
of the lu minosity density interior to the Su n. Like its pre- 
decessor llBinnev. Gerhard fc Soergell Il997l) . this model is 
based on the COBE L-band photometry, but it incorporates 
constraints from the slight difference in the mean apparent 
magnitude of red-clump stars on either side of the Galactic 
centre, and from tracers of spiral arms in the disk. Its bar is 
thi nner than that in the model of Bi nne y et a l. (1997). 

IBissantz. Englmaier fc Gerhardl feooah (hereafter 
BEG) determine the pattern speeds of the bar and the 
spiral structure, and upgrade the BG model to a mass 
model by simulating the flow of gas in potentials obtained 
by assigning spatially constant mass-to-light ratios to 
the BG model. Maps of the density of CO and HI in 
the longitude- velocity plane [(l,v) plane] are derived and 
compared with the corresponding observational plots. BEG 
find that they can reproduce the observed (l,v) plots by 
assigning to the bar a pattern speed that agrees with an 
independent determin ation from th e kinematics of the 
solar neighbourhood jDehnenl 1200(f) . while assigning a 
significantly lower pattern speed to the spiral structure. 
The spiral structure is four-armed and its amplitude in the 
mass density is larger by a factor ~ 1.5 than its amplitude 
in the near-infrared luminosity density. This standard 
model explains the entire rotation curve inside a galacto- 
centric radius of ~ 0.6i?o (where Ro is the galactocentric 
radius of the Sun) without any dark matter. However, 
at the solar radius Ro, the model yields a circular speed 
that is ~ 40kms _1 too low. To make up the balance, an 
axisymmetric quasi-isothermal dark halo with a large core 
radius (r c = 1.34i?o) and a small central density is added, 
and the mass-to-light ratio of the luminous component is 
lowered by less than 10% compared to the standard model. 

The mass model developed by BEG is strongly sup- 
ported by the results of surveys for microlensing events 
in fields towards the Galactic centre. Once the mass-to- 
light ratio of the BG model has been determined, a map 
of the microlensing optical depth for clump stars can be 
produced without any further assumptions, and BG present 
such a map. In units of 10 -6 , the EROS collaboration re- 
port an optical dep th r 6 = 0.94 ± 0.3 at (l,b) = (2.5°, -4°) 
iAfonso et al.l I2003T) while BG predicted tq — 1.2 at this 
location. The MACHO co llaboration report r 6 = 2.17±o;3g 
at (l,b) = (1.5°, -2.68°) JPopowski et al.ll2005l) . while BG 
predicted re = 2.4 here. Any transfer of matter from the stel- 
lar component to the dark halo would impair this excellent 
agreement. Still further evidence of the correctness of the 
Basel model has been provided by an JV-body model that 
self-consistently generates the mass den sity inferred by BEG 
(Biss antz. Debattista fc Gerhardl 120 04V This model, which 
has no free parameters, successfuly predicts the proper mo- 
tions of stars in various bulge fields. When a very plausi- 
ble model of the mass function of bulge and disk stars is 



adopted, it also reproduces the distribution determined by 
the MACHO collaboration of the duration of individual mi- 
crolensing events. Thus the Basel mass model, and its negli- 
gible contribution from CDM to the inner Galaxy, has been 
very thoroughly validated. 

The Basel group is now producing an upgraded ver- 
sion of the Basel model (Englmaier, private communication) 
and the calculations that follow are based on this new ver- 
sion. The model without dark matter (short-dashed curve of 
Fig. is 12% heavier than the luminous part of the model 
with a dark halo (dotted curve of Fig. . In the latter, the 
quasi-isothermal dark halo has a core radius r c = 10.7 kpc 
and an asymptotic velocity v x = 235kms _1 . It is assumed 
that Ro — 8 kpc and that v c (Ro) = 220kms _1 . The model 
with dark matter does not match the bumps in the rotation 
curve as quite as well as the standard model without dark 
matter does. 



3 THE MILKY WAY IN MOND 

The halo introduced by BEG is not cuspy like the halos pre- 
dicted by CDM cosmology (Diemand et al. 2004, and ref- 
erences therein). A cuspy halo could not have been added: 
indeed, what drives BEG to assign to baryons essentially all 
matter within the corotation radius of the bar, and to en- 
hance the amplitude of the spiral structure outside the bar, 
is the size of the non-circular velocities that are apparent 
in the observed (l,v) plots: the dark halo is assumed to be 
axisymmetric, so if much mass is shifted from the baryons to 
the halo, the non-axisymmetric component of the overall po- 
tential is weakened, and the non-circular velocities are pre- 
dicted to be too small. Given that the structure of the BEG 
halo is unexpected in the CDM theory, we now ask whether 
it can be eliminated by assuming that MOND is the correct 
theory of gravity. A full answer to this question requires 
exte nsive numerical work to solve the non-linear field equa- 
tion jBekenstein fc Milgromll984f) for the modified potential 
generated by the Galaxy, including contributions from the 
bar and spiral structure. However, the non-axisymmetry is 
important only at R < ^Ro, where the effects of MOND 
are small, which guarantees that the problem arising in the 
dark matter framework will not arise in MOND. Therefore, 
as in standard Newtonian theory, the first step is to study 
the circular-speed curve that follows from the axisymmetric 
component of the density distribution. 

3.1 The rotation curve 

In cylindrical or spherical symmetry the gravitational force 
per unit mass Kmond predicted by MOND in the Galactic 
plane is related to the corresponding Ne wtonian force p er 
unit mass KNowton by Milgrom's formula (Milgrom 1983): 

^(^mond/ao) Kmond = KNcwton, (1) 

where the interpolating function /i runs smoothly from 
(i(x) = i at i « 1 to (i(x) = 1 at i > 1. In a 
fiat axisymmetric disk, Milgrom's formula is onl y exact if 
MON D is viewed as a modification of inertia jMilgroml 
Il994f) rather than a modification of gravity. However, us- 
ing the field equation for the modified gravitational poten- 
tial, iBrada fc Milgroml jl995h have shown that equation Q 
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Figure 1. Circular-speed curves for the Milky Way. Dotted 
and short-dashed curves: Basel model with (dotted) and without 
(short-dashed) a quasi-isothermal dark halo. Long-dashed and full 
curves are MONDian curves for ao = 1.2 X 10 — 8 cm s -2 and either 
standard (short-dashed) or simple (full) interpolating functions. 
The dot-dashed curve is the MONDian circular-speed curve with 
the standard interpolating function and ao = 3.9 X 10 -8 cms -2 . 
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Figure 2. Terminal velocity curves in the fourth quadrant for 
the same models as in Fig. 1, compared with the data of Kerr 
et al. (1986). The terminal lign-of-sight velocities are defined by 
v t (l) = sign(l)v c (Ro sini) - v c (Ro)sinl. 



holds everywhere outside Kuzmin discs and disc-plus-bulge 
generalizations of them, while it provides a very good ap- 
proximation for exponential discs. As a consequence, this 
formula has been used to fit t he rotation curves of an imp res- 
sive list of external galaxies Zanders McOau ghll2002h . It 
is thus worthwhile to test the ability of this formula to fit 
the rotation curve of the Milky Way. 



Table 1. For the Basel model with and without dark mat- 
ter, and our 4 MONDian models, the columns display re- 
spectively the form of the interpolating function fi, the 
value of ao in units of 10 — 8 cm s — 2 , the value of the L-band 
mass-to-light ratio T^ in Mq / Lq (Tj, = £ 2 , where £ is the 
scaling constant for the rotation curve in BEG), the circu- 
lar speed at the solar radius Ro and at infinity in kms -1 , 
and finally the value fj,(xo) where xq = v 2 (Ro)/(Ro &o)- 
The values are displayed for models with Rq = 8 kpc. 
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Once the values of ao and of the mass-to-light ratio 
T are known, equation Q predicts the force field for each 
choice of interpolating function. From a sa mple of external 
galaxi es with high quality rotation curves iBeeeman et alJ 
<199lft derived a = 1.2±0.27x 1CT 8 cms -2 using the 'stan- 
dard' interpolating function fj,(x) = x/Vl + x 2 . Retaining 
the mass-to- light ratio of the Basel model without dark mat- 
ter, and using the standard interpolating function yields the 
long-dashed circular-speed curve shown on Fig. Q (Model I). 
The circular speed is then 207kms _1 at R = ^Ro, exactly 
as in the Basel model with dark halo, and 200 km s _1 at Ro- 
Thus making gravity MONDian with the standard interpo- 
lating function reduces the deficit in v c (Ro) from 40 kms 
to 20 kms" 1 , but does not eliminate it. 

However, the value of v c (Ro) is hard to determine 
and a lower value than 220 km s~ is not excluded (e.g. 
lOlling fc Merrifieldll99il . The primary observable is the ter- 
minal velocity vt(l) at each longitude /: 

v t (l) — sign(l)v c (Ro sin I) — v c (Ro) sini. (2) 

From vt(l) it is easy to determine v c (R) if one knows Ro 
and v c (Ro), but neither parameter is known with precision. 
Given this uncertainty, it is important to understand the 
predictions of each model for the run of v t (l) (see Fig. 0. 
When v c (Ro) = 220 kms -1 is assumed, as it was by the 
Basel group, the baryonic Basel model yields an excellent 
fit to the data at \l\ < 40°, but the curve is not self- 
consistent because it does not pass through zero at \l\ = 90°. 
The short-dashed curve in Fig. shows that when we self- 
consistently adopt v c (Ro) = 180 kms -1 (see Table 1), in the 
fourth quadrant the terminal velocities near / = —20° be- 
come marginally too negative and from there rise too steeply 
as / — > —90°. However, this self-consistent baryonic Newto- 
nian model is not clearly incompatible with the measured 
terminal velocities. 

Our MONDian Model I (long-dashed curve in Fig. 2) re- 
duces the discrepancy with the therminal velocity data, but 
still fits the data less well than the Basel model with dark 
matter (dotted curve in Fig. 2) because v c (Ro) = 200 kms - 
is low. If we enforce v c (Ro) = 220 kms" 1 simply by increas- 
ing T, we make v c too large at R ~ |i?o- A better fit can 
be obtained by adjusting both ao and T. To increase v c (Ro) 
we increase ao and thus trigger a MONDian correction at 
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smaller radii. This in turn obliges us to decrease T in order 
to keep a velocity of 207 km s -1 at R = ^Ro, exactly as in 
the Basel model with dark matter. This decrease in T re- 
quires a further increase in ao to retain v c (Ro) = 220 kms -1 , 
with the result that by the time we have achieved a satisfac- 
tory fit to the rotation curve at R = ^R and R = Ro, ao 
has increased significantly. The dot-dashed curves in Figs. 
and [5] shows the fits obtained when ao = 3.4 x 10 -8 cms -2 
and T L = 0.91 M Q /L Q (Model II). This model predicts an 
asymptotic circular speed of 210 kms - that is higher than 
those predicted by the conventional value of ao (see Table 
!)• 

The asy mptotic circular spe e d of t h e Galaxy is not well 
deter mined feinnev fc Dehnenl Il997t IWilkinson fcEvanj 
Il999l) . By contrast, in external galaxies with extended HI, 
the behaviour of v c at large R is strongly constrained 
by observations. Extended rotation curves have been ob- 
tained f or many galaxies with the aim of probing dark ha- 
los (e.g. iBegem an et alJll99ll : lBroeils fc van WoerderjlT99l 
iGentile et alj2004h . To be specific, we focus here on the case 
of NG C 3198, which was carefully studied bv lBegeman et all 
lll99llk this system is the textbook example of a galaxy with 
an extended flat rotation curve because it has a large appar- 
ent diameter and a velocity field that shows the gas to be 
confined to a plane and following accurately circular orbits. 
Fig. [3] shows that when one fits the rotation curve of NGC 
3198 with ao = 3.4 x 10 -8 cms -2 , the decrease in T that is 
required to match the measured asymptotic velocity results 
in v c being too low at small radii. 

Consequently, the data for NGC 3198 rule out the large 
value of ao required by Model II. In principle we can make 
Model II compatible with the standard value ao = 1.2 x 
10 -8 cms -2 by rescaling it: the rotation curves of the models 
listed in Table 1 are invariant if we change the value of Ro 
from 8 kpc and then multiply T by Ro/ (8 kpc) and ao by 
(8kpc)/i?o- However, to obtain ao = 1.2 x 10 -8 cms -2 by 
rescaling Ro, we would have to set Ro — 15 kpc, which is 
excluded by a wealth of data on the scale size of the Milky 
Way. 

From this we conclude that, with the standard inter- 
polating function, MOND cannot fit with perfect accuracy 
the rotation curves of both the Basel model a nd NGC 3198. 
This c onclusion is consistent with the work of IGentile et alJ 
i2004h . who concluded that by varying only T it is not possi- 
ble to obtain entirely satisfactory fits to the rotation curves 
of five spiral galaxies. 

We therefore consider an alternative 'simple' interpo- 
lating function, namely 



This function provides a less sudden transition from the 
Newtonian to the MONDian regime than does the stan- 
dard function. The continuous curve in Fig. shows that 
the simple interpolating function together with the conven- 
tional value ao = 1.2 x 10 -8 cms -2 yield a extremely good 
fit to the rotation curve of NGC 3198. 

The full curves in Figs. and [5] show the correspond- 
ing fits to the Basel model (Model III). We now have 
Vc(Ro) = 208 kms -1 . Although this value lies very close 
to the conventional value, we find it is impossible to push 
Vc(Ro) up to 220 kms -1 by increasing ao, because we then 
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Figure 3. Fits to the rotation curve of NGC 3198 using the pa- 
rameters of Model I (long-dashed), II (dot-dashed) and III (con- 
tinuous) of Table 1. Models I and III are compatible with NGC 
3198, while Model II is incompatible. This means that the value 
of ao used in Model II with the classical interpolating function 
is too high to account for a fit of NGC 3198. On the other hand 
the interpolating function of Model III provides an even better 
fit than the interpolating function used in Model II. 



have to decrease the mass-to-light ratio radically in order to 
fit the inner rotation curve [w c (0.5 Ro) = 207 kms -1 ], mak- 
ing Vc(Ro) too low again. The transition from Newtonian 
to MONDian physics provided by the simple interpolating 
function is insufficiently abrupt. Slightly changing the value 
of Ro does not eliminate the problem. However, Fig[5]shows 
that the fit of Model III to the terminal- velocity curve is ex- 
tremely satisfactory. This simple interpolating function thus 
fits the data for both NGC 3198 and the Milky Way. How- 
ever, it does not exactly reproduce the run of v c (R) of the 
Basel model; to achieve this, one needs a more complex form 
of the interpolating function. 



3.2 The interpolating function 

In low-surface-brightness (LSB) galaxies, x < 1 at all radii, 
with the consequence that the measured acceleration v^/R 
is approximately equal to [aoTGL(R)/ R 2 ] 1 ^ 2 . It follows that 
these galaxies constrain only the product aoT and leave 
the interpolating function unconstrained. In high-surface- 
brightness (HSB) galaxies such as NGC 3198, x ranges from 
values greater than unity down to small values. Hence with 
these objects we can break the degeneracy between ao and 
T, but, as Fig. |3] illustrates, the large gradients in rotation 
curves at small radii, combined with the poor spatial res- 
olution of obsevations in the 21-cm line, leave considerable 
degeneracy between ao and T, and constrain the interpo- 
lating function only weakly. The Basel model is strongly 
constrained at small radii, and it is pertinent to probe the 
extent to which it constrains the interpolating function in 
equation fT). 

Imagine an ideal galaxy in which both the asymptoti- 
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cally flat rotation curve and the luminosity distribution are 
known with precision. Then with x — v 2 / (Rao) we can write 



H(x) = T 



fjapx) 
clqx 



(4) 



where f(v%/R) is obtained from the Newtonian equation for 
the force per unit mass with the mass-density p(x) replaced 
by the luminosity-density j(x). From data at large x (corre- 
sponding to the central parts of the galaxy), where p(x) ~ 1 
and f(y) oc y, we can read off T. At large radii (small x) we 
have n{x) — x and f(y) oc y 2 , so we can determine do. Once 
T and ao are known, we can completely determine n(x) from 
equation @. 

As an example, we treat the circular-speed curve of the 
Basel model plus dark halo as a perfect set of data. From 
the centre we find T_l = 1.08 Mq/Lq, while the asymptotic 
velocity = 235 kms" 1 implies ao = 4.47 x 10 _8 cms~ 2 . 
This large value of ao suggests that the model's asymptotic 
velocity is too high. If one could show that the Milky Way's 
asymptotic circular speed really is this high, MOND would 
have been falsified. 

However, if we combine Tl — I.OSMq/Xq with the 
conventional value ao = 1.2 x 10~ 8 cms -2 , we infer an 
asymptotic velocity 170 kms -1 , and can read off the form of 
the interpolating function that perfectly fits the Basel model 
at R < Ro (Model IV). The dotted curve in Fig. 0] shows 
this function, which transitions smoothly from x/(l + x) at 
x ^ 1 to x/y/1 + x 2 at x > 10. 

Any successful underlying theory for MOND must pre- 
dict an interpolating functio n that lies betwee n the curves 
labelled I and III in Fig. H] iBekensteinl l)2004h introduced 
the TeVeS theory of gravitation, in which the physical met- 
ric depends on the Einstein metric, a vector field, and a 
scalar field. The dynamics of the scalar field is controlled by 
a Lagrangian density involving an undetermined function F 
that yields a MONDian behaviour for the potential in the 
nonrelativistic limit. Adopting a computationally tractable 
form for F, Bekenstein shows for the case of spherical sym- 
metry that Eq. (1) holds with p replaced by a function that 
depends on both the Newtonian acceleration and the gradi- 
ent of the scalar field. For accelerations that are not large 
compared to ao, the MONDian and Newtonian accelerations 
are then related by Eq. (1) with the interpolating function 



fi(x) = 



Vl + 4x- 1 



(5) 



V1 + 4Z + 1 

which is shown as the dashed curve in Fig.0] Since this curve 
lies far below the dotted curve at all values of x, the correc- 
tion to Newtonian gravity to which this function gives rise 
is much too big to be compatible with the dynamics of the 
Milky Way (see Table 1). Of course, since the Galaxy is not 
spherical, the interpolating function (5) does not apply to it, 
and even the concept of an interpolating function cannot be 
used for the Galaxy in TeVeS. However, it is not obvious that 
when non-spherical geometry is taken into account, Beken- 
stein's toy F will yield a sufficiently rapid transition from 
the MONDian to the Newtonian regime. Hence, groups that 
are endeavouring to determine the implications of TeVeS 
for large-scale structure should be cautious about conclu- 
sions obtained using Bekenstein's toy F, and especially not 
disregard this gravitational theory if predictions happen to 
disagree with observations when using this particular F. 




Figure 4. Fit (dotted line) to the interpolating function of 
MOND if one takes the Basel model as an ideal set of data from 
R = until Ro, a nd if one assumes an asymptotic velocity of 
170 kms -1 . The long-dashed and continuous lines correspond to 
the interpolating function of Models I and III, respectively. The 
short-dashed line is the interpolating function explored by Beken- 
stein (2004) 



3.3 The vertical equilibrium 

We now examine the consistency of MOND with the vertical 
dynamics of the solar neighbourhood. This is in principle a 
complex problem, because out of the plane it is not a priori 
clear that the Newtonian and modified-gravity accelerations 
will be parallel, as Milgrom's formula Q implies. However, 
equation Q is known to hold outside Kuzmin and exponen- 
tial disks, as well as in cylindrical symmetry, so it may be 
expected to provide a good first approximation in the Milky 
Way. 

At Ro (assumed here to be 8 kpc), the surface densit y 
of baryons is Eo = 53 Mq pc~ jHolmberg fc FlvnnlEooil . 
An infinite sheet with this surface density generates a New- 
tonian vertical acceleration -K^Nowton — 7rGEo — 0.25ao, for 
ao = 1.2 x 10 -8 cms -2 . The radial acceleration is v 2 /Ro — 
6.5 x K z Newton, which implies that the modification of the 
vertical force is dominated by the external field effect. Thus 
K z mond — -K Z Newton/£i(v c /-Ro ao). Observations show that 
1.1 kpc above the plane, the vertical force is larger than that 
arising from the baryons in Newtonian theory by a facto r 
~ 7/5 llKuiiken fc Gilmordll99ll iHolmberg fc FlvnrJl2004f) . 
It follows that MOND requires that 



lx{vl/Ro a ) ~ 5/7 ~ 0.71. 



(6) 



The extreme right column of Table 1 shows that, given 
the significant uncertainties, this constraint is satisfied by 
the standard, simple and fitted interpolating functions when 
a = 1.2 x 10~ 8 cms -2 . The choice a = 3.4 x 10~ 8 cms -2 
used in Model II places the solar neighbourhood so deeply 
into the MOND regime that the predicted large correction 
to the vertical force is excluded by the observations. 
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4 DISCUSSION 

By combining near-infrared photometry and simulations of 
gas flow in the plane, and without invoking dark matter in 
the inner 5 kpc, the Basel group has developed an extremely 
successful model of the Milky Way that accounts for the 
structure of (l,v) plots for CO and HI, for the proper mo- 
tions of bulge stars, for the microlensing optical depth to- 
wards bulge fields, and for the observed distribution of the 
durations of microlensing events. Given the number of these 
checks, there can be little doubt that we now really do know 
the distribution of baryons inside the solar radius. For no 
other galaxy do we have information of comparable quality. 

It is far from clear that the Basel model is compatible 
with the predictions of CDM. In the light of early indica- 
tions that baryons dominate the inner Galaxy, an attempt 
was made to build models that minimize the dark halo con- 
sistent with constraints from simu lations of clustering CDM 
jKlvpin. Zhao fc Somervillei2002h . In all these models, how- 
ever, CDM contributes significantly to the density at ~ 3 kpc 
from the Galactic centre where the Basel model requires the 
density to be almost entirely invested in stars. To investigate 
this matter further, we need a CDM-inspired model that in- 
cludes a stellar bar and reproduces the photometry of the 
Galaxy. It would be of great interest to calculate the predic- 
tions of such a model for microlensing surveys and the (I, v) 
diagrams of CO and HI. If those predictions disagree with 
obser vations, the CDM p aradigm would be strongly weak- 
ened. [Gentn^^d] i2004h have shown that the paradigm en- 
counters similar difficulties reproducing the rotation curves 
of a sample of five external galaxies. 

In MOND, there are two free parameters (T and oo, 
the latter a constant of nature) and one free function (//, 
to be fixed by an underlying theory of MOND). Once those 
quantities are known, the gravitational force field is com- 
pletely determined by the baryon distribution. We have in- 
vestigated what circular-speed curves MOND predicts from 
the distribution of baryons in the Basel model. A full an- 
swer to this question requires extensive numerical work, but 
a good approximation to the truth can be obtained from Mil- 
grom's formula Q . For the value of ao that has been deter- 
mined from observations of external galaxies, we encounter 
difficulty making v c (Rq) as large as the value, 220 kms -1 , 
assumed by the Basel model. When the standard interpo- 
lating function n(x) is used, satisfaction of this condition 
requires ao = 3.4 x 10 -8 cms -2 (Model II), a value which is 
incompatible with observations of NGC 3198. Keeping the 
conventional value ao = 1.2 x 10 _8 cms -2 (Model I) yields 
Vc(Ro) = 200 km s -1 implying slightly too low values for 
the terminal velocities. Interestingly, Gentile et al. (2004) 
have also found that MOND cannot account for the rota- 
tion curves of their sample of spiral galaxies, at least with 
that value of ao and that choice of interpolating function. 
However, in the Milky Way, a very good fit to the terminal 
velocity curve (Model III) can be obtained if the interpo- 
lating function (3) is used. It would be of great interest to 
test the ability of this interpolating function to fit the ro- 
tation curves of external galaxies in the sample of Gentile 
et al. (2004). We then showed how a perfect set of data 
could determine T, ao and the interpolating function n(x), 
and displayed the results (Model IV) that would follow if 



data were precisely those predicted by the Basel model at 
R ^ Ro. 

In the absence of dark matter, we conclude that the 
Galaxy's asymptotic velocity must be smaller than the 
value, 235kms _1 , implied by the dark halo of the Basel 
model (see Table 1); MOND predicts it to be ~ 170 ± 
5kms _1 . The fitted, numerically specified interpolating 
function of Model IV produces a local circular speed of 
220kms _1 , but lower values are obtained for the standard 
and simple interpolating functions of Models I and III, which 
suggests that 220 km s -1 is an upper limit for v c (Ro). We 
have also shown that it is unlikely that the interpolating 
function explored by Bekenstein (2004) as a toy model can 
account for the dynamics of the Milky Way. We finally 
showed that the three models with ao = 1.2 X 10 -8 cms -2 
(Models I, III and IV) are compatible with what is known 
about the vertical force at Ro — 8 kpc. 

It can thus be argued that, with a suitable choice of 
MOND's interpolating function, and a realistic choice of 
mass-to-light ratio, the Milky Way can be added to the long 
list of galaxies for which Milgrom's formula Q is success- 
ful in predicting the rotation curve from the baryon dis- 
tribution. The next step towards exploring the agreement 
between the data and the predictions MONDian gravity 
mak es from the Basel model i nvolves using a potential solver 
(e.g. Brada & Mikrom 1999) for the nonlinear Bekenstein- 
Milgrom equation to determine the MOND force generated 
by the non-axisymmetric components of the Basel model 
with the parameters and interpolating functions of our Mod- 
els I, III and IV. Another test for the relevance of MOND 
as an alternative to dark matter in the Milky Way will be 
provided by the measurement of the velocity dispersions in 
distant globular clusters orbiting the Galaxy in the outer 
halo, well into the deep-MOND regime (Baumgardt, Grebel 
& Kroupa 2005). 
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